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Abstract 

Background: The fungal pathogen Setosphaeria turcica causes turcicum or northern leaf blight disease on maize, 
sorghum and related grasses. A prevalent foliar disease found worldwide where the two host crops, maize and 
sorghum are grown. The aim of the present study was to find genes controlling the host defense response to this 
devastating plant pathogen. A cDNA-AFLP approach was taken to identify candidate sequences, which functions were 
further validated via virus induced gene silencing (VIGS), and real-time PCR analysis. Phylogenetic analysis was 
performed to address evolutionary events. 

Results: cDNA-AFLP analysis was run on susceptible and resistant sorghum and maize genotypes to identify 
resistance-related sequences. One CC-NB-LRR encoding gene GRMZM2G005347 was found among the up-regulated 
maize transcripts after fungal challenge. The new plant resistance gene was designated as St referring to 5. turcica. 
Genome sequence comparison revealed that the CC-NB-LRR encoding St genes are located on chromosome 2 in 
maize, and on chromosome 5 in sorghum. The six St sorghum genes reside in three pairs in one locus. When the 
sorghum St genes were silenced via VIGS, the resistance was clearly compromised, an observation that was 
supported by real-time PCR. Database searches and phylogenetic analysis suggest that the St genes have a 
common ancestor present before the grass subfamily split 50-70 million years ago. Today, 6 genes are present in 
sorghum, 9 in rice and foxtail millet, respectively, 3 in maize and 4 in Brachypodium distachyon. The St gene 
homologs have all highly conserved sequences, and commonly reside as gene pairs in the grass genomes. 

Conclusions: Resistance genes to 5. turcica, with a CC-NB-LRR protein domain architecture, have been found in 
maize and sorghum. VIGS analysis revealed their importance in the surveillance to S. turcica in sorghum. The St 
genes are highly conserved in sorghum, rice, foxtail millet, maize and Brachypodium, suggesting an essential 
evolutionary function. 



Background 

The immune system has developed in a stepwise manner 
by progressive sophistication of basic functions that helped 
ancestral organisms to survive in their hostile environ- 
ment. Recognition of pathogens in a species-specific way 
results in the generation of a very robust mode of surveil- 
lance system in plants. This form of protection termed 
resistance (R) protein-mediated or effector-triggered 
immunity is induced when a plant encoded R protein 
"perceives" the presence of a pathogen-derived effector 
molecule, represented by specific avirulence (Avr) gene 
products [1]. Following recognition of the pathogen, one 
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or more signal transduction pathways are induced in the 
host plant and these lead to the prevention of colonization 
by the pathogen. 

The majority of characterized R proteins encode a 
nucleotide-binding site (NB) and leucine-rich repeats 
(LRR). NB-LRR-encoding genes make up one of the lar- 
gest and most variable gene families found in plants, with 
most plant genomes containing several hundred family 
members [2-6]. The N-terminal ends of R-proteins are 
predominantly composed of a TIR (Toll/Interleukin-1 
Receptor) homologous domain or form a coiled-coil (CC) 
motif. Monocots in particular, have numerous CC-NB- 
LRR proteins in their genomes. Accumulating data suggest 
furthermore that N termini of R-proteins may interact 
with a range of pathogen-derived proteins. However, the 
LRR domain may determine the final outcome of this 
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recognition, leading to downstream signaling and initiation 
of defense responses [7] . 

Many R-genes are located in clusters that either com- 
prise several copies of homologous sequences arising from 
a single gene family or co-localized R-gene sequences 
derived from unrelated gene families [8,9]. This genomic 
make-up assists multiple proteins to become modified via 
various genie and intergenic processes enabling rapid evo- 
lution and adaptation to changes in a pathogen genome 
[10]. R-genes can also act in pairs [11,12]. The R-gene 
pairs can differ in genomic location and protein domain 
structure but also to their interaction with different patho- 
gen isolates. 

The heterothallic ascomycete Setosphaeria turcica (Lut- 
trell) Leonard & Suggs (anamorph: Exserohlium turcicum, 
former Helminthosporium turcicum) causes turcicum or 
northern leaf blight disease on maize. This fungal patho- 
gen also attacks sorghum and related grass species, for 
example Johnson grass [13,14]. Turcicum leaf blight is one 
of the most prevalent foliar diseases in most maize-grow- 
ing regions of the world. The disease causes periodic epi- 
demics associated with significant yield losses, particularly 
under conditions of moderate temperature and high 
humidity [15-17]. Resistance to S. turcica has mainly been 
characterized in maize. S. turcica was earlier named 
Helminthosporium turcicum and resistance has hitherto 
been designated Ht and conferred by major race-specific 
genes {Htl, Ht2, Ht3 or HtN) or via partial resistance, 
reviewed by Welz and Geiger [18]. In our work we 
designate the new resistance genes as St referring to 
Setosphaeria turcica. 

Maize and sorghum are the most important staple cer- 
eals for sub-Saharan Africa (SSA). While maize is an 
introduced crop [19], sorghum is believed to have been 
domesticated in SSA particularly in the Nile basin or 
Ethiopia, as recently as 1000 BC [20] . Sorghum like many 
other crop species experience large problems with plant 
pathogens, particularly fungal diseases. Turcicum leaf 
blight incited by S. turcica is one main problem [21]. 
This disease has been considered as of minor importance 
in Uganda until 1988 when it caused extensive yield 
losses on maize [22]. By introducing improved resistance 
in new varieties the threat posed by the disease was sub- 
sequently reduced. Severe and sporadic outbreaks of tur- 
cicum leaf blight have now reappeared in East Africa 
[23-25]. A change in the S. turcica population has been 
suggested to be the main cause of this shift in disease 
pattern. In order to detect potential new changes of the 
S. turcica pathogen and the turcicum leaf blight disease, 
a survey was undertaken in Uganda to examine the sor- 
ghum - S. turcica pathosystem in terms of disease sever- 
ity and incidence, race patterns and new resistant 
resources [26]. It can be concluded from those studies 
that fungal isolates from sorghum could infect maize. 



Upon cross inoculation on maize differential lines har- 
boring different Ht genes, four S. turcica isolates were 
identified as race 1, two as race 2, and one isolate corre- 
sponded to race 0 and race 3, respectively, whereas 10 
isolates were unclassified. Highly resistant sorghum 
accessions originating from a regional collection were 
also identified. 

In this work, we used cDNA-amplified fragment length 
polymorphism (AFLP) on resistant and susceptible maize 
and sorghum genotypes to identify differentially expressed 
genes, when challenged with S. turcica. This was followed 
by functional assessment of selected gene candidates by 
virus-induced gene silencing (VIGS) using a Brome mosaic 
virus vector. We found one R-gene cluster, containing six 
CC-NB-LRR encoding genes residing as three pairs in the 
sorghum genome, of importance for defense to S. turcica. 
Genome data further showed that the St genes are highly 
conserved within monocots. 

Results 

Identification of an up-regulated /?-gene family in maize 
and sorghum in response to S. turcica inoculation 

In order to identify important defense genes to S. turcica, 
cDNA-AFLP analysis was carried out on susceptible (S) 
and resistant (R) sorghum and maize genotypes following 
fungal infection. In our case, the Ugandan sorghum gen- 
otypes GA06/18 (R) and Sila (S) and the maize A619Htl 
(R) and A619 (S) lines were used. The sorghum material 
had earlier been evaluated on various agronomical traits 
including important fungal diseases. Apart from S. tur- 
cica responses, GA06/18 was found to be susceptible to 
Cercospora sorghi, and Colletotrichum sublineolum, 
whereas Sila was susceptible to C. sorghi and resistant to 
C. sublineolum. 

In total, approximately 3000 transcript-derived frag- 
ments were monitored ranging from 50 to 600 bp in size 
using different primer combinations (Additional file 1). 
Unique, up- or down-regulated transcripts in the resistant 
genotypes compared to the susceptible, sampled at 24, 48 
and 72 hours post inoculation (hpi) were excised, ampli- 
fied, sequenced and analyzed for putative function. The 
final transcript-set comprised of 68 sorghum and 82 maize 
gene candidates. Among these genes, 11 and 13, respec- 
tively, were putative stress-related according to closest 
genes identified in other organisms using BLASTP. 

One CC-NB-LRR encoding putative R-gene 
(GRMZM2G005347), a member of a homologous gene 
pair with GRMZM2G005452 in the same locus on chro- 
mosome 2, and uniquely expressed in the resistant maize 
genotype, was further studied (Figure ID). Genome analy- 
sis revealed presence of 6 homologous genes in sorghum 
(Figure 1A). These six genes were given the prefix St refer- 
ring to S. turcica and designated StlA (Sb05g008280), StlB 
(Sb05g008140), St2A (Sb05g008350), St2B (Sb05g008030), 
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Figure 1 Evolution of St R-gene cluster in monocots species. Chromosome location, duplication and ancestry of the St gene cluster in A. 
sorghum, B. millet, C. rice, D. maize, and E. Brachypodium. Events preceding, (light grey) and post speciation (dark grey) are shown. F. Proposed 
ancestral R-gene cluster composition using an ancestral tree of grass species adapted from Bowman et al. [55], and phylogenetic analysis of 
homologous genes in each species. Genes colored in relation to St genes as follows; Stl blue, St2 red, and St3 green. Gene information is listed 
in Additional file 2. 



St3A (Sb05g0082S0), and St3B (Sb05g008270). Quantita- 
tive real-time PCR confirmed furthermore that five (StIA, 
St2A, St2B, St3A and St3B) of the six St genes showed 
high relative transcript levels when the sorghum resistant 
GA06/18 plants were challenged with S. turcica (Figure 2). 
One gene, StlB, was expressed to a much lower extent 
compared to the other St genes, outside the limit of detec- 
tion. In Sila, only St2B and St3A showed a significant 
increase (P < 0.005) in expressions when challenged with 
5. turcica (Figure 2). 

The St genes are conserved among grasses 

The six St genes in sorghum form three gene pairs in a 
cluster on chromosome 5 and share a common ancestor 
(Figure 1; Additional file 2; Additional file 3). St gene 
orthologs were also found in clusters when searching 
the rice, maize, foxtail millet and Brachypodium genome 
databases. The St gene encoded proteins from the other 
grass species, grouped with the sorghum St proteins 
with high edge support (100) (Additional file 2). The 
rice genome contains orthologs of sorghum StIA, StlB, 
St2A f St2B and an St3 gene (Figure 1A, B). This indi- 
cates that the ancestor of rice and sorghum likely had a 
copy of these genes. Sorghum St3A and St3B are likely a 
result of a more recent genome duplication event after 
the split between the rice and sorghum species (Figure 
1G). The rice genome also contains multiple copies of 




GA06/18 



Figure 2 Relative qPCR values of StIA, St2A, St2B, St3A and 
St3B transcripts in sorghum GA06/18 and Sila plants 
inoculated with S. turcica, 24 hpi. Water treatment was used on 
respective genotypes as a control. Error bars indicate standard 
deviation between three biological replications. * Indicates a 
significant increase {P < 0.05) compared to control levels. Primer 
used listed in Additional file 2. 
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StlA, St2A and St2B orthologs, likely produced from 
gene duplications after the species split from sorghum. 
Likewise, the Setaria italica (foxtail millet) genome con- 
tains orthologs of StlA, StlB, St2A and St2B, with seven 
of the nine genes found in a cluster within the same 
scaffold, as complete chromosome annotation have yet 
to be determined (Figure 1C). An St3 homolog was not 
found in millet. In addition to the maize gene pair iden- 
tified in our cDNA-AFLP analysis, BLASTP and 
BLASTN searches revealed a third single gene homolog, 
GRMZM2G050959, St2A on maize chromosome 2 (Fig- 
ure ID). The model grass Brachypodium genome, on 
the other hand, has a gene pair orthologous to StlB on 
chromosome 4, and one to St2B on chromosome 5, but 
lacks all other gene homologs (Figure IE). The St gene 
cluster is maintained between sorghum, rice and possi- 
bly millet genomes but is smaller in maize and Brachy- 
podium with St genes located across or on different 
chromosomes. 

Sequence homology was also found between sorghum St 
proteins and Arabidopsis CC-NB-LRR encoding genes 
(Figure 3; Additional file 4). All six St proteins formed a 
cluster together with the CC rather than TIR domain con- 
taining R proteins from Arabidopsis indicating a closer 
evolutionary relationship as expected. The nearest related 
Arabidopsis gene is RPM1, a gene mediating resistance to 
Pseudomonas syringae isolates expressing the avrRpml or 
avrB genes [27]. 




Figure 3 Un-rooted maximum-likelihood phylogram inferred 
from nucleotide binding (NB) and leucine rich repeat (LRR) 
domains, of six resistance proteins to S. turcica (St) in 
sorghum, compared with Arabidopsis NB-LRR resistance 
proteins with known function. LR-ELW values above 75% are 
shown. A. NB-LRR resistance proteins with a coiled-coil (CC) domain 
at the N-terminal end. B. NB-LRR resistance proteins with a Toll/ 
lnterleukin-1 receptor (TIR) at the N-terminal end. Units indicate 
substitutions/site. R-proteins used are listed in Additional file 4. 

V J 



Adapting the VIGS system on sorghum 

Genetic transformation of sorghum and maize is possible 
but laborious and requires other genotypes than those 
used in this study to be successful [28,29]. Hence, our 
candidate genes were further studied using virus induced 
gene silencing (VIGS) using the Brome mosaic virus 
(BMV) system, previously used to silence genes in mono- 
cots [30]. VIGS was followed by fungal inoculation to 
assess the potential defense function of the St genes. In 
our hands, the VIGS procedure was not successful when 
applied to the A619Htl maize genotype. Because the St 
genes were up-regulated upon fungal inoculation with 
S. turcica in our sorghum GA06/18 genotype (Figure 2), 
we continued the studies on our sorghum materials. 

Two VIGS constructs (1 and 2) with high identity to 
the 6 St genes in sorghum were designed (Figure 4) 
including examination for their off-target gene silencing 
capacity. The highest non-St sorghum gene similarity 
belongs to a related i?-gene pair, Sbl0g028720 and 
Sbl0g028730, located in a different subgroup upon phy- 
logenetic analysis (Additional file 2), and used as a con- 
trol for off-target gene silencing. The selected sequences 
were amplified and ligated into the third plasmid 
(pF13m) in the BMV system, and used to infect the sor- 
ghum plants. 
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Figure 4 Schematic alignment of VIGS constructs and St genes 
and their closest off-target genes. Two sections of StlA were PCR 
amplified and ligated into VIGS plasmid pF13m. The identity of the 
constructs to each of the St genes is shown. E values for construct 
1 and 2 against St and closest off target genes; StlA 2e-128, 2e-128; 
St IB 4e-68, NA; St2A 2e-34 5e-48; St2B 1e-43 6e-60; St3A 2e-33 8e- 
27; St3B 3e-32 4e-24; Sbl0g028730 3e-07, NA; Sbl0g028720 3e-06, 
NA. NA indicates no significant similarity. 
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The VIGS procedure was first optimized. Sorghum 
seeds were surface sterilized before sowing to minimize 
additional stress by other microorganisms. mRNA was 
produced by in vitro transcription, added to inoculation 
buffer and rubbed directly onto the second leaf of three 
week old sorghum plants. No intermediate step involving 
barley as virus host was used. The virus spreads systemi- 
cally throughout the plant with silencing greatest in the 
second and third leaves above the inoculation site and 
complete silencing rarely achieved [30]. Seven days post 
infection (dpi), light green colored streaks were visible on 
the third leaf, indicating viral symptoms and successful 
infection by the virus. In order to confirm onset of silen- 
cing quantitative real time-PCR was carried out on leaf 
samples from the VIGS treated plants (Figure 5). There 
was a significant decrease in the relative transcript levels 
in relation to control plants inoculated with empty plas- 
mid suggesting a clear down-regulation of five of the six 
targeted genes, particularly by construct 1, in both sor- 
ghum genotypes. Relative transcript levels of 
Sbl0g028720 and Sbl0g028730 were not influenced in 
VIGS treatments indicating no off- target silencing. 

Silencing of St genes increases S. turcica infection in the 
resistant and susceptible sorghum genotypes 

Fungal colonization and growth on plants inoculated with 
the different VIGS constructs compared with control 
material was carefully monitored. The different phenotypic 
observations are summarized in Figure 6; and Additional 
file 5. Fungal growth was further assessed by detaching 
infected leaves and placing them in a petri dish containing 
moist filter paper followed by incubation in the dark at 
25°C for two days, as described by Levy [31]. The develop- 
ment of conidiophores protruding through leaf lesions fol- 
lowed by rapid asexual spore development indicated 
fungal colonization of the leaf material, and a susceptible 
phenotype. 

A hypersensitive response (small dark/red spots) 
occurred at 2 dpi on the resistant GA06/18 genotype 
upon fungal challenge while the plants treated with 
empty vector produced a somewhat delayed HR pheno- 
type 3 dpi. When VIGS construct 1 was applied to 
GA06/18 plants prior to fungal inoculation, larger and 
more numerous lesions with chlorotic halos developed 
compared to the control plants. Disease lesions spread 
laterally along the leaf and fungal conidiophores and 
spores were produced under sporulating conditions. 
Similarly, when the effect of construct 2 was assayed, 
the disease lesions were seen 2 dpi and spread laterally 
to form large lesions that produced large numbers of 
fungal spores. The disease lesions were larger than those 
induced by construct 1, at 7 dpi. On the susceptible Sila 
plants clear disease symptoms, necrotic spots, and 
chlorotic halos around fungal appressoria were seen 
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Figure 5 Relative qPCR values of St gene transcripts when 
inoculated with VIGS constructs 1 and 2 compared to empty 
vector control. Both GA06/18 (A) and Sila (B) genotypes showed 
down-regulation of five of the six St genes when inoculated with 
either one, or both of the constructs. Error bars indicate standard 
deviation between three biological replications. * Indicates a 
significant decrease {P < 0.05) compared to control levels. 



2 dpi. Large numbers of asexual fungal spores were pro- 
duced on conidiophores protruding from necrotic 
lesions. When Sila plants were inoculated with the 
empty VIGS vector, prior to fungal inoculation, similar 
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Figure 6 Leaf phenotypes of resistant GA06/18 and susceptible Sila plants, 7 dpi with S. turcica, pretreated either with water, empty 
vector or VIGS construct 1 and 2 



disease symptoms occurred 2 dpi. In contrast, on Sila 
plants inoculated with our VIGS construct 1, slightly 
larger and more frequent lesions appeared compared to 
control plants. The disease symptoms were further 
amplified when construct 2 was used, resulting in larger 
necrotic lesions, and profuse fungal sporulation. In 
order to correlate these observed disease phenotypes 
with fungal growth, fungal DNA was quantified in the 
VIGS materials (Figure 7). S. turcica DNA increased to 
1.5 ± 0.4 pg/ng sorghum DNA in GA06/18 leaves 
inoculated with VIGS construct 1, and to 3.6 ± 0.9 pg/ 
ng sorghum DNA when using construct 2, from a near 
zero level in control plants (non-VIGS and empty vec- 
tor). A significant (P < 0.005) increase in fungal DNA 
was also found in samples from Sila inoculated with 
construct 1 (1.2 ± 0.4 pg/ng sorghum DNA), and con- 
struct 2 (0.8 ± 0.9 pg/ng sorghum DNA), compared to 
control samples with approximately 0.5 pg/ng sorghum 
DNA. 

Taken together, as expected the resistant GA06/18 gen- 
otype showed a compromised defense response when 
inoculated with VIGS construct 1 or 2 prior fungal inocu- 
lation. Interestingly, we observed enhanced disease pheno- 
types on the susceptible Sila genotype upon corresponding 
VIGS treatments. 



Discussion 

Sorghum [Sorghum bicolor (L.) Moench] serves as a 
major food staple and fodder resource especially in arid 
and semi-arid regions of the world [32]. It is mainly a 
self-pollinating and diploid grass species (2n = 2x = 20), 
with a genome size of 1C = 730 Mbp, which is about 25% 
the size of the maize genome [4,5]. In the sorghum gen- 
ome, 211 NB-LRR encoding 7?-genes are present, which 
is approximately half the number found in rice and 
slightly more compared to Arabidopsis [4]. The number 
of NB-LRR encoding genes in the small genome of the 
wild grass Brachypodium is estimated to 178 [6]. But in 
the much larger maize genome, 95 NB-LRR encoding 
genes have up to now been identified [33]. However, 
depending on search programs and threshold settings, 
slightly different 7?-gene numbers in each grass species 
are published. 

It is postulated that the high numbers of i?-genes in 
plant genomes and their large sequence diversity are 
essential evolutionary factors in the surveillance machin- 
ery to resist pathogen attacks. Resistance genes evolve 
through duplication, unequal crossing over, recombina- 
tion and diversification leading to clusters of paralogous 
genes [10,34]. The proliferation of i?-genes is also 
coupled with rapid turnover of gene copies, eventually 
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Figure 7 Real-time qPCR assessment of S. turcica DNA on 
resistant GA06/18 and susceptible Sila leaf samples, 7 dpi with 
S. turcica, pretreated either with water, empty vector, or VIGS 
construct 1 and 2. Error bars indicate standard deviation between 
three biological replications. * Indicates a significant increase {P < 
0.05) in fungal DNA compared to control levels. 
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leading to deletion or expansion and thus dynamic R- 
gene clusters [33]. Resistance gene clusters have also 
been found to be conserved between different species in 
Poaceae [35], although, such clusters are in the minority 
with 71.6% being specific to a species [33]. 

Whole genome duplications occurred when the grass 
subfamilies diverged from each other and genome data 
suggest further, that paleo-duplicated gene pairs in sor- 
ghum and rice remained extant in about 17% of the cases 
[36]. Recent duplications of chromosomal segments are 
particularly found on rice chromosomes 11 and 12, and 
corresponding regions on chromosome 5 and 8 in sor- 
ghum. Chromosome 5, in the sequenced BTx623 sor- 
ghum genotype, where the St genes are located showed 
the highest abundance (62) of i?-genes [4]. Thirty-six of 
these NB-LRR encoding genes are affected by recent 
duplication events based on the bioinformatic analyses 
presented by Wang et al. [36], including St3A and St3B, 
which is in agreement with our results (Figure 1; Addi- 
tional file 2). Interestingly, the rice genome contains 
orthologs of StlA, StlB, St2A, St2B and a single ortholog 
of the St3 genes, all in one single locus. This indicates 
that this gene cluster predates the species split of rice 
and sorghum. In the grass family, sorghum, maize and 
millets belong to the same sub-family (Panicoideae), 
whereas rice is located in Ehrhartoideae [37]. It is 



estimated that these two subfamilies diverged from a 
common ancestor 50-70 million years ago together with 
Pooideae, the subfamily to which Brachypodium, wheat, 
and barley belong. 

In a genome-wide comparison of Arabidopsis thaliana 
and A. lyrata, the evolutionary pattern of the i?-genes 
could be divided into two distinct groups, the positively 
selected (> 50%) with high sequence divergence between 
the two species, or the stably selected genes (< 30%) [38]. 
The remaining genes were only found in one genome and 
absent from the other. The St genes found in this work 
have experienced few sequence exchanges resulting in low 
divergence, and hence more resemble the description of 
stably selected genes, although the copy numbers vary 
between the five grass genomes compared (Figure 1). That 
NB-LRR encoded i?-genes remain conserved between dif- 
ferent grass species is presently believed to be a common 
phenomenon [33]. 

Sorghum plants, particularly genotypes with red seed 
color, accumulate a range of phenolic substances in 
response to pathogen attacks [39] . Large amounts of red- 
pigmented flavonoids induced at the site of infection were 
also seen in our materials, particularly in the resistant gen- 
otype. Whether flavonoids contribute to the defense 
response against S. turcica is not elucidated but a genetic 
link has been found in the sorghum - C. sublineolum 
interaction, produced via the presence of 3-deoxyantho- 
cyanidins [40]. Reinforcement of plant cells via callose 
deposition upon pathogen attacks have been observed in 
many pathosystems. Enhanced callose deposition has also 
been reported as a resistance response to S. turcica in 
maize [41]. Despite extensive staining efforts, no callose 
accumulation was seen in either of our sorghum genotypes 
(data not shown). 

Furthermore, our gene silencing work resulted in an 
enhanced susceptible response in Sila, our susceptible sor- 
ghum cultivar. This observation may suggest that by tar- 
geting the St genes in this genomic background, effects on 
downstream signaling masked in the resistant sorghum 
genotype are revealed, and could potentially constitute a 
fraction of the quantitative traits earlier found [41]. This 
hypothesis is speculative and remains to be included in 
future functional studies of the St genes. Future studies do 
also comprise a search for important effectors in the gen- 
ome recently released from JGI http://www.jgi.doe.gov. In 
parallel, the sequence information from the St gene cluster 
is presently converted into molecular markers and used in 
germplasm assessments and breeding programs in East 
Africa, an important development to sustain sorghum and 
maize crop production in this part of the world. 

Conclusions 

Our cDNA-AFLP analysis on susceptible and resistance 
maize and sorghum genotypes challenged by S. turcica 
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resulted in identification of a CC-NB-LRR encoding 
gene in maize. This gene resides in two loci on maize 
chromosome 2. In sorghum, 6 St orthologous genes are 
present in a cluster of three pairs, on chromosome 5. 
Upon gene-silencing of the sorghum St genes, the resis- 
tance was clearly compromised, an observation that was 
supported by real-time PCR analysis and fungal DNA 
quantification. Database searches and phylogenetic ana- 
lysis suggest that the St genes have a common ancestor 
present before the subfamily split, 50-70 million years 
ago, and the genes are highly conserved in sorghum, 
rice, foxtail millet, maize and Brachypodium. 

Methods 

Plant and fungal materials 

Resistant (R) and susceptible (S) Sorghum bicolor geno- 
types from Uganda, GA06/18 (R) and Sila (S), and maize 
lines A619Htl (R) and A619 (S) provided by USDA ARS, 
were used in the study. The plants were grown in a growth 
chamber (Percival) using a 12/12 h photoperiod at 22°C. A 
single spore isolate from S. turcica infected sorghum (Igl), 
or infected maize (Mbl), collected from Iganga and 
Mbale, Uganda, was used for all sorghum and maize analy- 
sis, respectively. The fungal DNA was extracted using a 
modified CTAB method [42] . DNA was analyzed by using 
S. turcica specific ITS1 and ITS2 primers (F -GCAA- 
CAGTGCTCTGCTGAAA and R- AT A AG ACG G CC A A- 
CACCAAG). PCR was carried out using the following 
conditions: 10 ng of template DNA was added to a 24 \i\ 
mix consisting of H 2 0, 2.5 mM MgCl 2 , 2.5 [i\ Taq buffer 
(Fermentas, Helsingborg, Sweden) 0.2 mM of each dNTP, 
0.25 (iM of forward and reverse primers and 1 U of Taq 
polymerase (Fermentas) with: 3 min at 94°C, 35 cycles of 
(1 min at 94°C, 1 min at 60°C, and 1.5 min at 72°C), and 
final extension at 72°C for 10 min. The PCR products 
were separated on 1% agarose gels to confirm fragment 
size, (344 bp) followed by sequencing (Macrogen Inc., 
Seoul, Korea). 

Fungal inoculation of plant material 

Three-week old seedlings were inoculated on the third leaf 
whorl with 25 \i\ conidia suspension (5 x 10 5 conidia/ml) 
as described by Carson [43]. Inoculated leaves from three 
to four plants were pooled and harvested at 24, 48, and 
78 hours post inoculation (hpi) for cDNA-AFLP analysis. 
Water treated control samples were harvested at the same 
time-points. 

RNA extraction and cDNA-AFLP analysis 

Total RNA was isolated from the leaf samples using the 
BioRad RNA isolation kit (BioRad, California, USA) fol- 
lowed by mRNA preparation with the mRNA capture kit 
(Roche, California, USA). cDNA was synthesized with 
Oligo-dT primer and RevertAid™ H Minus M-MuLV 



Reverse Transcriptase (Fermentas). Second strand was 
synthesized using E. coli DNA Polymerase I (Fermentas). 
The double stranded cDNA was digested with BstYl and 
Msel (Fermentas) and ligated to respective adaptors, pre- 
amplified and later selectively amplified using the BstYl 
+N ( 33 P labeled) and Msel +N primers. Pre-amplification 
was carried out with the adapter-ligated cDNA, Taq 
DNA Polymerase (Fermentas) and the non-selective pri- 
mers specific to the BstYl and Msel adapters using 25 
cycles of 94°C for 30 s; 56°C for 1 min and 72°C for 
1 min. The pre-amplified reaction mixture was diluted 
600-fold and 5 \i\ was used for final selective amplifica- 
tion with 24 primer combinations, carried out with BstYl 
+N ( 33 P labeled) primers (Additional file 1) and touch- 
down amplification [44]. The selective amplification pro- 
ducts were resolved on 6% polyacrimide gel run at 
100 W until 4300 Vh was reached. Gels were dried and 
exposed to Kodak Biomax film (Amersham Pharmacia, 
California, USA) for 5-7 days. 

Isolation and sequencing of transcripts 

Approximately 150 transcripts (unique, up and down- 
regulated) from the resistant genotypes in relation to the 
susceptible genotypes, were excised from the dried PAGE 
gels, eluted in H 2 0 and PCR amplified using the non- 
selective primers under the same conditions as earlier 
described in the pre-amplification step. The products were 
cloned into the pJET 1.2 blunt vector (CloneJET™ PC, 
Fermentas) and sequenced. The sequences were analyzed 
using the BLASTN and BLASTX programs [45] and com- 
pared with sequences deposited in NCBI, GRAMENE and 
PHYTOZOME databases. Identified fungal sequences 
were excluded. 

Virus induced gene silencing (VIGS) in sorghum 

The VIGS system used is based on the monocot-infecting 
Brome mosaic virus (BMV) as previously described [30] 
but pre-inoculation on barley was excluded. The BMV 
VIGS vector consists of three plasmids harboring BMV 
RNA1 (pl-1), RNA2 (p2-2) and RNA3 (pF13m, formally 
pF3-5/13 A/G ), respectively. To generate VIGS constructs, 
PCR fragments ranging from 246 to 253 bp in size were 
amplified from the sorghum candidate gene using geno- 
mic DNA of the resistant GA06/18 genotype and gene- 
specific primers harboring Ncol and Avrll restriction 
sites using the Primer 3 version 0.4.0 http://frodo.wi.mit. 
edu/primer3/ software (Additional file 6). Prior to PCR 
amplification, off-target gene searches were undertaken 
to design optimal VIGS constructs (Figure 4). After 
restriction, each fragment was cloned into the corre- 
sponding site of the pF13m plasmid. The identity of the 
inserts was verified by sequencing. Pl-1, p2-2 and the 
pF13m containing different constructs were digested 
with Spel, PshAl and Psh Al, respectively. Infectious RNA 
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transcripts were synthesized from linearized plasmids 
through in vitro transcription using T3 Polymerase (Fer- 
mentas), according to manufacturer instructions. 1 [A of 
the reaction product was run on a 1.5% agarose gel to 
confirm presence of a transcript. 

Plant inoculation procedures were performed as 
described [30] with slight modifications. A 10 \A aliquot of 
the transcription mix from each of the plasmids pl-1, p2- 
2 and pF13m-insert was combined with 30 (il FES inocula- 
tion buffer and used directly to rub inoculate the second 
and third leaves of 3 -week-old sorghum and maize plants. 
As a control, plants were inoculated in the same way with 
water or combined transcripts from pl-1, p2-2 and empty 
pF13m. Maize and sorghum plants were challenged with 
S. turcica as earlier described one week after viral inocula- 
tion (when faint chlorosis and vein clearing started to 
appear) to assess the effect of the different constructs. 
Plants were randomized and coded to reduce potential 
bias in the scoring of fungal colonization and growth. 

Quantitative real-time PCR 

Prior to fungal inoculation of the VIGS treated sorghum 
plants, approximated 3 cm of the second leaf above the 
VIGS inoculated leaf was collected from 3 independent 
plants in triplicates for each condition and used for RNA 
extraction as previously described. First strand-cDNA was 
synthesized from 1 (ig of total RNA, with Oligo-dT primer 
and RevertAid™ H Minus M-MuLV Reverse Transcrip- 
tase (Fermentas) according to the manufacturer's instruc- 
tions. Real-time PCR was carried out using the first strand 
cDNA in an iQ5 cycler (Bio-Rad). Maxima Sybr Green/ 
Fluorescein qPCR Master Mix (Fermentas) was used for 
PCR amplification in a 20 \A total reaction volume consist- 
ing of 10 ul of SYBR Green qPCR Master Mix, 0.3 uM for- 
ward and reverse primers and 5 ng of cDNA template. All 
PCRs were performed in triplicate under the following 
amplification conditions; 10 min at 95°C followed by 
40 cycles of 95°C, for 15 s, 60°C for 30 s, and 72°C for 
30 s, followed 1 min at 95°C, and melt curve analysis. Pri- 
mers sequences for St genes were designed using the Pri- 
mer 3 version 0.4.0 http://frodo.wi.mit.edu/primer3/ 
software (Additional file 7). The sorghum elongation fac- 
tor 1-alpha (Sb02g036420) and Actin (Sb01g010030) were 
used as reference genes and relative transcript values were 
calculated. All calculations and statistical analyses were 
performed as described in the ABI PRISM 7700 Sequence 
Detection System User Bulletin #2 (Applied Biosystems, 
USA) slightly modified by Vetukuri et al. [46]. Quantifica- 
tion of S. turcica DNA on VIGS material 7 days post fun- 
gal inoculation was carried out as earlier described [41]. 
Approximately 3 cm of leaf material from three plants was 
pooled and DNA extracted. Three biological samples per 
treatment were analyzed. Statistical significance was calcu- 
lated using Student's t-test. 



Genome analysis 

The amino acid sequences of SUA (Sb05g008280), StlB 
(Sb05g008140), St2A (Sb05g008350), St2B (Sb05g008030), 
St3A (Sb05g008250) and St3B (Sb05g008270) were aligned 
to sorghum, maize, millet, rice, Brachypodium and Arabi- 
dopsis genome databases using BLASTN and BLASTP 
(PHYTOZOME). Predicted domains were identified using 
coiled-coil prediction [47]), LRRfinder [48] and CD-Search 
[49]. St-hke gene loci were identified using Genomic 
Viewer (PHYTOZOME). Phylogenetic analysis was con- 
ducted using Treefinder and maximum likelihood and 10, 
000 replicates [50]. The JTT+G model [51] was found to 
best fit the data using ProtTest v2.4 [52]. Confidences 
were calculated using local rearrangement of expected 
likelihood weights (LR-ELW) [53]. Phylograms were 
drawn using Treeview 1.6.6 [54]. 

Additional material 



Additional file 1: PCR primer combinations used in cDNA-AFLP 
analysis. 

Additional file 2: Maximum likelihood phylogenetic tree using the 
model JTT+G based on amino acid sequence from the coiled coil 
(CC), nucleotide binding (NB) and leucine rich repeat (LRR) domains 
of St proteins in sorghum, and closely related R proteins. Names 
refer to PH^TOSOME gene identifier. Physcomitrella patens /?-protein 
Pp1s1_327V6, was used as an out-group. LR-ELW edge support values 
are shown [53], Substitutions per site are indicated. 

Additional file 3: Information on genes used in Figure 1 and their 
putative function. Data retrieved from the PHYTOSOME database. 
GenBank accession numbers are stated where present. 

Additional file 4: Information on Arabidopsis genes used in Figure 3. 
Data retrieved from the TAIR database. 

Additional file 5: Disease phenotypes on sorghum leaves, 
monitored 1-12 days post inoculation (dpi) with S. turcica on the 
resistant wild type GA06/18 and the susceptible Sila cultivar. The 

plants were treated with either water, empty BMV vector, construct 1 or 
construct 2, prior to fungal inoculation. The data is compiled from 25-30 
plants per BMV construct and controls. The experiment was repeated 2 
times. 

Additional file 6: Gene specific primers for VIGS constructs. 
Restriction sites are in bold. 

Additional file 7: List of primers used in real time PCR analysis. 



Acknowledgements 

We greatly acknowledge B. Sarosh for assisting with the cDNA-AFLP work. 
This study was funded by the SIDA-SAREC grant for the East African 
Biotechnology Network (BIO-EARN), Sida SWE-2005-453, and the Swedish 
University of Agricultural Sciences. 

Author details 

1 SLU, Uppsala Biocenter, Dept. Plant Biology and Forest Genetics, Uppsala P. 
0. Box 7080, S-750 07, Uppsala, Sweden. 2 Dept. of Crop Science, Makerere 
University P.O. Box 7062, Kampala, Uganda. 3 National Agriculture Research 
Organisation, P.O. Box 295, Entebbe, Uganda. 

Authors' contributions 

TM carried out the qPCR, genomic analysis, created figures, and ran the 
cDNA-AFLP and VIGS analyses together with MB. IF supported on the VIGS 
analysis. CD and PO conceived the study and participated in writing 



Martin et al. BMC Plant Biology 201 1, 11:151 
http://www.biomedcentral.eom/1471-2229/1 1/151 



Page 1 0 of 1 1 



together with all authors. All authors read and approved the final 
manuscript. 

Received: 3 June 201 1 Accepted: 3 November 201 1 
Published: 3 November 201 1 

References 

1. Jones JDG, Dangl J: The plant immune system. Nature 2006, 
444(71 17):323-329. 

2. Meyers BC, Kozik A, Griego A, Kuang H, Michel more RW: Genome-wide 
analysis of NB-LRR-encoding genes in Arabidopsis. Plant Cell 2003, 

15(4):809-834. 

3. Liu J, Liu X, Dai L, Wang G: Recent progress in elucidating the structure, 
function and evolution of disease resistance genes in plants. Journal of 
Genetics & Genomics 2007, 34(9)765-776. 

4. Paterson AH, Bowers JE, Bruggmann R, et al: The Sorghum bicolor genome 
and the diversification of grasses. Nature 2009, 457(7229):55 1-556. 

5. Schnable PS, Ware D, Fulton RS, et al: The B73 maize genome: complexity, 
diversity and dynamics. Science 2009, 326(5956):1 112-1115. 

6. The International Brachypodium Initiative: Genome sequence and analysis 
of the model grass Brachypodium distachyon. Nature 2010, 
463(7282):763-768. 

7. Collier SM, Moffett P: NB-LRRs work a "bait and switch" on pathogens. 
Trends in Plant Science 2009, 1 4(1 0):521 -529. 

8. Smith SM, Pryor AJ, Hurlbert SH: Allelic and haplotypic diversity at the rp1 
rust resistance locus of maize. Genetics 2004, 167(4):1 939-1 949. 

9. Kuang H, Woo S-S, Meyers BC, Nevo E, Michelmore RW: Multiple genetic 
processes result in heterogeneous rates of evolution within the major 
cluster disease resistance genes in lettuce. Plant Cell 2004, 
16(11):2870-2894. 

10. McDowell JM, Simon SA: Recent insights into R gene evolution. Molecular 
Plant Patholology 2006, 7(5)437-448. 

11. Dixon MS, Jones DA, Keddie JS, Thomas CM, Harrison K, Jones JDG: The 
tomato Cf-2 disease resistance locus comprises two functional genes 
encoding leucine-rich repeat proteins. Cell 1996, 84(3)451-459. 

12. Eitas TK, Dangl JL: NB-LRR proteins: pairs, pieces, perception, partners 
and pathways. Current Opinion in Plant Biology 2010, 13(4)472-477. 

13. Hamid A, Aragaki M: Inheritance of pathogenicity in Setosphaeria turcica. 
Phytopathology 1975, 65(3):280-283. 

14. Chiang M, van Dyke C, Leonard K: Evaluation of endemic fungi for 
potential biological control of Johansongrass {Sorghum halepense): 
Screening and host range tests. Plant Disease 1989, 73(6)459-464. 

15. Perkins JM, Pedersen WL: Disease development and yield losses 
associated with northern leaf blight on corn. Plant Disease 1 987, 
71(10):940-943. 

16. Carson ML, van Dyke C: Effect of light and temperature on expression of 
partial resistance of maize to Exserohilum turcicum. Plant Disease 1994, 

78(4)408-411. 

17. Pratt R, Gordon S: Breeding for resistance to maize foliar pathogens. Plant 
Breeding Reviews 2006, 27:1 19-173. 

18. Welz G, Geiger H: Genes for resistance to northern corn leaf blight in 
diverse maize populations. Plant Breeding 2000, 119(1):1-14. 

19. Mangelsdorf PC: Corn: its origin, evolution and improvement. Belknap 
Press, Cambridge, Mass. USA; 1974. 

20. Kimber CT: Origins of domesticated sorghum and its early diffusion into 
India and China. In Sorghum: origin, history technology and production. 
Edited by: Smith CW, Frederiksen RA. John Wiley 2000:3-98. 

21. Ngugi HK, King SB, Holt J, Julian AM: Simultaneous temporal progress of 
sorghum anthracnose and leaf blight in crop mixtures with disparate 
patterns. Phytopathology 2001, 91(8):720-729. 

22. Adipala E, Lipps EP, Madden LV: Occurrence of Exserohilum turcicum on 
maize in Uganda. Plant Disease 1993, 77(1 ):202-205. 

23. Ebiyau J, Oryokot OE: Sorghum {Sorghum bicolor (L.) Moench. Agriculture 
in Uganda. Crops: National Agricultural Research Organisation Fountain 
Publ; 200111. 

24. Muiru WM: Histological studies and characterization of races of 
Exserohilum turcicum the causal maize agent of northern leaf blight of 
maize in Kenya. PhD thesis University of Nairobi, Kenya; 2008. 

25. Ramathani I: Characterisation of turcicum leaf blight epidemics and 
pathogen populations in the Exserohilum turcicum - Sorghum 
pathosystem in Uganda. MSc thesis Makerere Univ. Kampala, Uganda; 2009. 



26. Ramathani I, Biruma M, Martin T, Dixelius C, Okori P: Disease severity, 
incidence and races of Setosphaeria turcica on sorghum in Uganda. 
European Journal of Plant Pathology 201 1, 131(3)383-392. 

27. Grant MR, Godiardt L, Straube S, Ashfield T, Lewald J, Satler A, Ines RW, 
Dangl JL: Structure of the RPM1 gene enabling dual specificity disease 
resistance. Science 1995, 269(5225):843-846. 

28. Opabode J: Agrobacterium-mediated transformation of plants: emerging 
factors that influence efficiency. Biotechnology and Molecular Biology 
Reviews 2006, 1(1): 12-20. 

29. Gurel S, Gurel E, Kaur R, Wong J, Meng L, Tan HQ, Lemaux PG: Efficient, 
reproducible 4grobacfe/7um-mediated transformation of sorghum using 
heat treatment of immature embryos. Plant Cell Reports 2009, 
28(3)429-444. 

30. Ding XS, Schneider WL, Chaluvadi SR, Mian MAR, Nelson RS: 
Characterization of a Brome mosaic virus strain and its use as vector for 
gene silencing in monocotyledonous hosts. Molecular Plant-Microbe 
Interactions 2006, 19(1 1):1 229-1 239. 

31. Levy Y: Variation of fitness among field isolates of Exserohilum turcicum 
in Israel. Plant Disease 1991, 75(2):1 63-1 66. 

32. Doggett H: Sorghum. John Wiley, New York, US;, 2 1988. 

33. Li J, Ding J, Zhang Y, Wang J, Chen J-Q, Tian D, Yang S: Unique 
evolutionary pattern of numbers of gramineous NBS-LRR genes. 
Molelucar Genetics & Genomics 2010, 283(5)427-438. 

34. Meyers BC, Kaushik S, Nandety RS: Evolving disease resistance genes. 
Current Opinion in Plant Biology 2005, 8(2):1 29-1 34. 

35. Luo S, Peng J, Li K, Wang M, Kuang H: Contrasting Evolutionary patterns 
of the Rp1 resistance gene family in different species of Poaceae. 
Molecular Biology & Evolution 201 1, 28(1)313-325. 

36. Wang X, Wang H, Paterson AH: Seventy million years of concerted 
evolution of a homoeologous chromosome pair, in parallel, in major 
Poaceae lineages. Plant Cell 201 1, 23(1):27-37. 

37. Bolot S, Abrouk M, Masood-Quraishi U, Stein N, Messing J, Feuillet C, Salse J: 
The inner circle of the cereal genomes. Current Opinion in Plant Biology 
2009, 12(2):1 19-125. 

38. Chen Q, Han Z, Jiang H, Tian D, Yang S: Strong positive selection drives 
rapid diversification of R-genes in Arabidopsis relatives. Journal of 
Molecular Evolution 2010, 70(2):1 37-148. 

39. Dicko MH, Gruppen H, Barro C, Traore AS, van Berkel WJH, Voragn AGJ: 
Impact of phenolic compounds and related enzymes in sorghum 
varieties for resistance and susceptibility to biotic and abiotc stresses. 
Journal of Chemical Ecology 2005, 31(1 1):2671 -2687. 

40. Ibraheem F, Gaffoor I, Chopra S: Phytoalexin-dependent resistance to 
anthracnose leaf blight requires a functional yellow seed! in Sorghum 
bicolor. Genetics 2010, 184(4):91 5-926. 

41. Chung C-L, Longfellow JM, Walsh EK, Kerdieh Z, van Esbroeck G, Balint- 
Kurti P, Nelson RJ: Resistance loci affecting distinct stages of fungal 
pathogenesis: use of introgression lines for QTL mapping and 
characterization in the ma\ze-Setosphaeria turcica pathosystem. BMC 
Plant Biology 2010, 10:103. 

42. Okori P, Rubaihayo PR, Ekwamu A, Fahleson J, Dixelius C: Genetic 
characterization of Cercospora sorghi from cultivated and wild sorghum 
and its relationship to other Cercospora fungi. Phytopathology 2004, 
94(7):743-750. 

43. Caisson ML: Inheritance of latent period length in maize infected with 
Exserohilum turcicum. Plant Disease 1995, 79(6):581-585. 

44. Vos P, Hogers R, Bleeker M, Reijans M, Vandelee T, Homes M, Frijters A, 
Pot J, Peleman J, Kuiper M, Zabeau M: AFLP a new technique for DNA 
fingerprinting. Nucleic Acids Research 1995, 23(21)4407-4414 

45. Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment 
search tool. Journal of Molecular Biology 1990, 215(3)403-410. 

46. Vetukuri RR, Avrova AO, Grenville-Briggs LJ, van West P, Soderbom F, 
Savenkov El, Whisson SE, Dixelius C: Evidence for involvement of Dicer- 
like, Argonaute, and Histone Deacetylase proteins in gene silencing in 
Phytophthora infestans. Molecular Plant Pathology 201 1, 12(8):772-785. 

47. Lupas A, Van Dyke M, Stock J: Predicting coiled coils from protein 
sequences. Science 1991, 252(501 0):1 162-1 164. 

48. Offord V, Coffey TJ, Werling D: LRRfinder: a web application for the 
identification of leucine-rich repeats and an integrative Toll-like receptor 
database. Development & Comparative Immunology 2010, 34(1 0):1 035-1 041. 

49. Marchler-Bauer A, Bryant SH: CD-Search: protein domain annotations on 
the fly. Nucleic Acids Research 2004, 32A/V327-331. 



Martin et al. BMC Plant Biology 201 1, 11:151 
http://www.biomedcentral.eom/1471-2229/1 1/151 



Page 11 of 1 1 



50. Jobb G, von Haeseler A, Strimmer K: TREEFINDER: a powerful graphical 
analysis environment for molecular phylogenetics. BMC Evolutionary 
Biology 2004, 4:18. 

51. Jones DT, Taylor WR, Thornton JM: The rapid generation of mutation data 
matrices from protein sequences. Computer Applications in the Biosciences 
1992, 8(3):275-282. 

52. Abascal F, Zardoya R, Posada D: ProtTest: selection of best-fit models of 
protein evolution. Bioinformatics 2005, 21 (9):21 04-21 05. 

53. Strimmer K, Rambaut A: Inferring confidence sets of possibly misspecified 
gene trees. Proceedings of the Royal Society B: Biological Sciences 2002, 
269(1 487):1 37-1 42. 

54. Page RD: TreeView: an application to display phylogenetic trees on 
personal computers. Computer Applications in Biosciences 1996, 
12(4):357-358. 

55. Bowman J, Floud S, Sakakibara K: Green genes - comparative genomics of 
the green branch of life. Cell 2007, 129(2):229-234. 



doi:1 0.1 1 86/1 471 -2229-1 1-151 

Cite this article as: Martin et al.: A highly conserved NB-LRR encoding 
gene cluster effective against Setosphaeria turcica in sorghum. BMC Plant 
Biology 2011 11:151. 



Submit your next manuscript to BioMed Central 
and take full advantage of: 

• Convenient online submission 

• Thorough peer review 

• No space constraints or color figure charges 

• Immediate publication on acceptance 

• Inclusion in PubMed, CAS, Scopus and Google Scholar 

• Research which is freely available for redistribution 

Submit your manuscript at (^\ RioMM i r pnt ^\ 

www.biomedcentral.com/submit Blomea ^ enirai 



